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^ ! Abstract 



Carrying out a research program outlined by John S. Bell in 1987, we arrive 
at a relativistic version of the Ghirardi-Rimini-Weber (GRW) model of spou- 



ts , taneous wavefunction collapse. The GRW model was proposed as a solution of 

^Q I the measurement problem of quantum mechanics and involves a stochastic and 

vQ ' nonlinear modification of the Schrodinger equation. It deviates very little from 

(^ . the Schrodinger equation for microscopic systems but efficiently suppresses, for 

^^ I macroscopic systems, superpositions of macroscopically different states. As sug- 

gested by Bell, we take the primitive ontology, or local beables, of our model to be 
a discrete set of space-time points, at which the collapses are centered. This set 
is random with distribution determined by the initial wavefunction. Our model is 
P3 ' nonlocal and violates Bell's inequality though it does not make use of a preferred 

5 ' slicing of space-time or any other sort of synchronization of spacelike separated 

qh| points. Like the GRW model, it reproduces the quantum probabilities in all cases 

presently testable, though it entails deviations from the quantum formalism that 
are in principle testable. Our model works in Minkowski space-time as well as in 
(well-behaved) curved background space-times. 



PACS numbers: 03.65.Ta; 03.65.Ud; 03.30. +p. Key words: spontaneous wave- 
function collapse; relativity; quantum theory without observers. 



1 Introduction 

. . . / am particularly struck by the fact that the [GRW] model is as Lorentz 
invariant as it could be in the nonrelativistic version. It takes away the 
ground of my fear that any exact formulation of quantum mechanics must 
conflict with fundamental Lorentz invariance. J. S. Bell |S] 
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In 1986, Ghirardi, Rimini, and Weber (GRW) proposed a model of spontaneous wave- 
function collapse ^ni based on a stochastic and nonlinear modification of the Schrodinger 
equation. When combined with a clear ontology, the GRW model turns quantum me- 
chanics into a completely coherent theory. It resolves all paradoxes of quantum mechan- 
ics, in particular the measurement problem, and accounts for all phenomena of quantum 
mechanics in terms of an objective reality governed by mathematical laws. However, 
the GRW model is nonrelativistic. 

In 1^ and again in jHj, John S. Bell emphasized that the GRW model has a property — 
multi-time translation invariance — that can be regarded as a nonrelativistic surrogate 
of Lorentz invariance. This fact suggests that the biggest difficulty one would expect 
with turning a nonrelativistic theory into a relativistic one — the difficulty caused by the 
lack of a temporal ordering of spacelike separated events — is absent in the GRW model 
right from the start. We find this suggestion to be correct, and can indeed specify a 
relativistic version of the GRW model. We proceed along the lines of Bell's suggestions; 
in particular, we do not use a continuous spontaneous collapse model (corresponding to 
a diffusion process in Hilbert space and known as continuous spontaneous localization, 
or CSL), but rather a discrete one corresponding to a jump process in Hilbert space. 
Furthermore, we follow Bell in taking as the primitive ontology, or local beables, of the 
model the space-time points where the collapses are centered jH El E|. "A piece of 
matter then is a galaxy of such events" j3|. We will call these points "flashes." 

For a recent overview of spontaneous collapse models, see [3j. We regard our model as 
a step towards one possible explanation of the probability rules of quantum theory in the 
relativistic realm. This realm differs in two ways from that of nonrelativistic quantum 
mechanics, to which the GRW model applies: in the requirement of Lorentz symmetry, 
and in the phenomenon of particle creation and annihilation typical of quantum fleld 
theory. Here we will not be concerned with the latter, but focus on covariance, and thus 
on relativistic quantum mechanics; correspondingly, what we shall mean by "relativistic" 
is "Lorentz invariant," or its analogue in curved space-time. 

Our relativistic model is surprisingly similar to the original GRW model, which it 
approaches in the nonrelativistic limit. Its structure is in no way more complicated than 
that of the GRW model in Bell's flash-based version. The two models have the following 
features in common: (i) the only objects in the universe (beyond the given space-time 
geometry) are the wavefunction and the flashes; (ii) two new constants of nature are 
needed, the collapse rate 1/r per particle and the width a of the localization; (iii) time 
reversal invariance is broken, while (in flat space-time) rotation, space translation, time 
translation, parity, and gauge invariance are obeyed; (iv) the dynamics is intrinsically 
stochastic. 

Our model is based on relativistic quantum mechanics of A^ particles. The question of 
identical particles we plan to address in a separate work; here we shall avoid this question 
and base our considerations on the quantum mechanics of distinguishable particles. We 
shall use the letter i to denote (apart from \/—T) the particle types, i G {!,..., A^}, and 
Qi for the set of flashes belonging to i; the elements of this set are timelike separated 
from each other. 



The wavefunction is a multi-time wavefunction, i.e., it is defined on the Cartesian 
product of A^ copies of space-time. We use the Dirac equation as the relativistic ver- 
sion of the Schrodinger equation determining the evolution of the wavefunction apart 
from the collapses (but we will mostly not worry whether the wavefunction lies in the 
positive energy subspace, except in Section ITTj) . More precisely, we use the multi-time 
formalism with A^ Dirac equations. For the consistency of this set of equations, we 
cannot have interaction potentials. To avoid discussing the question of interaction in 
relativistic quantum mechanics, we will assume non-interacting particles. Interaction 
can presumably be included by allowing for particle creation and annihilation, which 
however is beyond the scope of this paper. In any case, the difficulty of including inter- 
action that we encounter here does not stem from the spontaneous collapses but rather 
from the mathematics of multi-time equations, and is thus encountered by every kind 
of relativistic quantum mechanics. 

We will give three equivalent descriptions of our model. In one of them, we refer to 
an arbitrary slicing (foliation) of space-time into spacelike surfaces and obtain a Markov 
process for the temporal evolution (relative to this slicing) of wavefunction and flashes; 
since this picture is not manifestly covariant, we postpone it to the end. We will begin 
instead with an iterative construction of the flashes, and take this to be the deflnition 
of our model. In another description, we provide a formula for the joint distribution of 
the flashes on space-time in terms of the initial wavefunction. 

This paper is organized as follows. In Section |21 we recall the deflnition of the GRW 
model. In Section 13 we deflne our relativistic variant. In Section |3 we compute joint 
distributions of the flashes and flash rates. In Section El we reformulate the model in 
terms of a temporal evolution relative to an arbitrary spacelike slicing of space-time. 
In Section IHl we show that the low velocity limit of our model is the GRW model. In 
Section [7| we discuss some predictions of our model. Finally, in Section |H1 we conclude 
by comparing it to other models in the literature. 

2 The GRW Model 

2.1 Definition 

We briefly recall the GRW model, following Bell's description 4J. The wavefunction 
v]/ = \|/(ri, . . . , rjv, t) evolves unitarily between the collapses. At the time T when a 
flash of type I G {1, . . . , A^} occurs, at location X G M^, the wavefunction collapses 
according to 



1/2 



pr{x,T- 



where the jump factor is a Gaussian 



i(r)=irexp(-^) (2) 



whose width a is a new constant of nature, of order of magnitude 10 ^ni, and the 
normahzation constant K is chosen so that 



d'r\3{r)\' = l. (3) 

Furthermore, 



-') 



Pi{x,t)= d^ri---d^rN\j{ri-x)^{ri,...,rN,t)\\ (4) 

so that the collapsed wavefunction in (^ is normalized. The rate for a collapse of type 
i to occur in the volume element d^y is 

-Pi{y)d^y, (5) 

r 

where r is another new constant of nature, of order of magnitude 10^^ sec. To put this 
differently, 

FToh(Qin[t,t + dt]xR' = {Y},Y e[t,t + dt]xd'y\^t) = ^^^{^t\Uyf\^t), (6) 

where ji{y) is the self-adjoint "collapse" operator that multiplies by the function j{ri — 
y), and Qi is the set of all flashes with label i (a subset of space-time). 

2.2 Another Definition 

We now give a second, equivalent, formulation of the GRW model that is closer to the 
way we will formulate the definition of our relativistic model. A first change consists in 
that, rather than having the wavefunction change discontinuously at some time T, we 
will speak of two different wavef unctions, one representing the situation before collapse 
and the other the situation after (or better, one uncollapsed and one collapsed), and 
extend both wavefunctions to all times, future and past, using the unitary, collapse-free 
evolution. Hence, in this terminology and notation, wavefunctions never collapse. If 
that seems paradoxical, we emphasize that the goal is to define the distribution of the 
random sets Qi, as it is only the flashes that chairs, tables, and observers are made of, 
and whatever mathematical formulation leads to the right distribution is allowable. 

Assume, as we will do in the relativistic model, that the A^ particles do not interact, 
that is, that the Hamiltonian is of the form H = Hi + . . . + Hj^j where H^ acts only 
on the i-th coordinate of the wavefunction — so that Hi commutes with Hj, i ^ j. This 
allows us to define the Schrodinger evolution also for a multi-time wavefunction 

^(ri, ti, . . . , r^, tM) = e-^^^*^/" ■ ■ ■ e-^^-*-/^^(ri, 0,...,r^,0). (7) 

The building block of the reformulation of the GRW model is the following procedure 
for obtaining, from given time values Ti, . . . , T/v (later taken to be times of flashes) and 



a given wavefunction \1/ on R^^, new random time values T[, . . . , T^ (later taken to be 
times of subsequent flashes), associated random locations Yi, . . . , Yjy G M^, and a new 
(collapsed) wavefunction $ on R^^: Let ATi, . . . , AT^r be independent, exponentially 
distributed random variables with expectation r, and set T/ = Tj + ATj. The joint 
distribution of the Y^, 



ti 



|2 



N)\ ■ 

(9) 



Prob(ri e d^y^, . . . , l"7v e £yj^) = p{y^, ■■■,yN) d^Vi ■ ■ ■ d^VN, 
has density p : M^^ -^ M defined by 

piyi,---,yN) = / d^Zi---d^ZN \j{zi - y^) ■ ■ ■ j{zN - yN)'^izi,T{, . . . , zn,T'^ 

Now define, for all (zi, . . . , z^) G M^^, 

$(zi,T„...,z^,T^)- pi/2(l-,,...,r^) ^10) 

and extend $ to all other times values by unitary multi-time evolution such as ((Tj). 

The random sets Qi of flashes are obtained by iterating this procedure. We start 
with Ti = 0, . . . , Tjv = and the initial wavefunction, obtain the first flash (T/, Yi) in 
every Qi, and then take the times of these flashes and the new wavefunction as the input 
of the next round of the procedure. It can be shown that this definition is equivalent to 
the definition given in the previous subsection as it leads to the same random sets Qi. 
We also obtain in this way a sequence of wave functions, each on M^^, each a solution 
to (|7j), and each associated with a set of A^ flashes. 

3 The Relativistic Model 

In this section we define the relativistic model. We will use a rather abstract formulation 
that has the advantages of being concise, simple, manifestly consistent, and manifestly 
covariant. 

3.1 The Dirac Equation 

We begin with recalling some relevant aspects of the Dirac equation and introducing 
some notation along the way. We generally use the letters x,y, . . . ,X,Y, . . . to denote 
space-time points, where capital letters usually stand for random space-time points. 

As the particles are not interacting, our model does not require that they live in the 
same space-time; instead, we may have N space-time manifolds ^i, . . . , ^n (which may 
be flat or curved), such that the wavefunction is a function on ^i x ■ ■ ■ x ^^ = Yli ^i, 
and the set Qi of flashes of type z is a discrete subset of ^i. We find that this greater 
mathematical generality, though physically unnecessary, facilitates the mathematical 
treatment of the model. 



The evolution of the wavefunction \1/ = "^{xi, . . . , x^) is determined, apart from the 
collapses, by the Dirac equation 

^H{Vi,^ - fA,^{x,))^ = m,^. (11) 

Here, rrii and Cj are mass and charge of particle i, Vj is the (covariant) derivative on 
^i, and Ai the electromagnetic vector potential on ^j. \E' takes values in (C^)®^ or, 
in the case of curved space-times, is a cross-section of the vector bundle 

^= U %i,...,x^) = U ^l,xi®---®^7V,x^, (12) 

where ^j is the bundle of Dirac spin-spaces associated with ^j. 

As in Section \'2.'2\ we will, throughout Section 13 take the wavefunctions to always 
obey the collapse-free evolution, this time given by the multi-time Dirac equation (jllll . 
and to be defined on all of Yli -^i- In particular, at any collapse we will consider two 
wavefunctions rather than one that changes discontinuously. 

A remark is necessary on consistency of the multi-time formalism. Multi-time equa- 
tions are not always consistent; what can go wrong is that the propagator for the i-th 
time coordinate (relative to some coordinate system) fails to commute with the j-th 
propagator, j ^ i, in which case it is impossible to find, for every initial wavefunction, 
a wavefunction on Yl- ^i solving all N evolution equations simultaneously. The set of 
equations ()lip. however, is consistent because the metric and vector potential acting on 
the coordinate Xi do not depend on the other coordinates Xj, j 7^ i. 

We now briefly explain how the Dirac equation (jllj) defines unitary propagators on 
suitable Hilbert spaces. Consider first the one-particle Dirac equation 

z^7^(V^-fA^)^ = m^. (13) 

We need and assume from now on that there are no closed timelike curves. With every 
spacelike surface S is associated a Hilbert space I>^(S), and with any two spacelike 
Cauchy^ surfaces S and S' is associated a unitary propagator JJ^ : -^^(S) -^ L'^iX'') 
as follows. i^^(S) contains cross-sections of the bundle & of Dirac spin spaces (which 
on Minkowski space-time are just C^) restricted to S and is endowed with the scalar 
product 

{m\^) = I d^x^{x)n^,{x)-i^^{x) (14) 

where n^(x) is the (future-directed) unit normal vector on S at x, and the volume 
measure d?x is the one arising from the Riemann metric on S. For simplicity, we will 
often write 

|^(,T)|2=^(x)n^(a;)7^^(x) (15) 

when it is clear which surface S we are considering. Unitarity of the propagator U^ 
follows from the continuity equation 



VJ^7'^* =0, (16) 



"'^A Cauchy surface is a surface intersected by every complete timelike curve. 



which is a consequence of flT^ . Of course, U^, U^ = U^ , and tj^ is the identity op- 
erator on L^(S). For multi-time wavefunctions, we will only consider spacelike surfaces 
of product form, Si x ■ ■ ■ x Stv- Then, the associated Hilbert space coincides with the 
tensor product of the Hilbert spaces associated with the Sj, and the unitary propagator 
between two such spaces coincides with the tensor product of the unitary propagators 
for each i. 

3.2 Notation 

We need some more notation. Let ^{x) denote the future of x, i.e., the future light 
cone and its interior including x itself, and <^(x) the past of x. When S* is a subset of 
space-time we write ^{S) ("the future of S"') for Uxe5^('^) ^^^ ^{S) ("the past of 
5"') for |Jxe5 ^{x). For y G ^(x), let t-dist(?/,a;) be the timelike distance of y from x, 
i.e., the supremum of the lengths of all timelike curves connecting x to y; for Minkowski 
space-time, 

t-dist(y,x) = {{y^ - x^){y, - x,)y'\ (17) 

For r > 0, let 

^(x) = {ye ^{x) : t-dist(?/, x) = r} (18) 

be the future surface of timelike distance r from x, or the r-hyperboloid; for r = this 
is the future light cone. For y G <^{x), we also write Jif{y,x) = Jift-dist{y,x){x) for the 
hyperboloid centered at x containing y. If S is a spacelike surface and x, y G S, we 
write s-dists(x, y) for the spacelike distance from x to y along S, i.e., the infimum of the 
Riemannian lengths of all curves in S connecting x to y. The speed of light is denoted 
by c. 

3.3 Law of the Flashes 

We are now prepared to write down the stochastic law of the flashes. The building block 
is the following procedure for obtaining, from given flashes Xi, . . . , Xj\f and a (normal- 
ized) wavefunction \E^ obeying (lllj). new flashes Yi, . . . , Yj^j and a new wavefunction $ 
obeying (fTT|. 

Let ATi, . . . , ATjsf be independent, exponentially distributed random variables with 
expectation r. Choose {Yi, . . . ,Yn) at random from Si x ■ ■ ■ x S^r = Ilj^*' where 
Sj = J^cATii^i), with distribution 

Prob(ri G d^yi, . . . , F^ G d^yN) = pivi, ■■■^Vn) d^Vi ■ ■ ■ d^VN (19) 

as follows. The volume of d^yi is computed using the Riemann metric on Sj; the distri- 
bution density p : flj Sj — ;> M is deflned by 



p{yi,---,yN)= I dZi---dZN\JT.^{yi^Zi)---j's:j^{yN,ZN)'^{zi,...,ZN)\] (20) 



and, for any spacelike surface S, the jump factor js : S x S ^ M is defined by 

,,fe.^) = if.Wexp(-!:^5!||^) (21) 



with Kj](z) chosen so that 



Now define $ on Yli ^i by 



d'y\My,z)\' = l. (22) 



*^"^' • • • ' ^^) pV2(r,,...,y^) ' (23) 

and define it on the remainder of Yl^ ^i by extending via (jllll . 

To the extent that the hyperboloids S, = ^at,(^j) are Cauchy surfaces (see remarks 
in the subsequent subsection), $ is uniquely determined on Yli^i from initial data 
on Yli^i (and normalized due to (UHl)), ^ is normalized on Hj^i) ^^d, by (j221), P is 
normalized in the sense 

/ d^yi---d^yNp{yi,...,yN) = l. (24) 

Note also that F^ G ^(X^). 

Our relativistic model is defined by iterating this procedure. As initial data, specify 
a wavefunction \|/ = v|/0 and one flash Xf of every type i. Apply the above procedure 
to obtain \l/^ = $ and a new flash X^ = YJ of every type. Repeat the procedure with 
\1' = \1'^ and Xi = X}, and so on. In this way, obtain for every type i a random sequence 
of flashes, Qi = {X^, X}, Xf, . . .}. 

Concerning the initial flashes of the universe, it seems an idea worth considering that 
the Big Bang, i.e., the initial singularity of the space-time geometry, is the space-time 
location of the initial flash for each i G {1, . . . , N}. This could have the status of a law, 
and would remove the arbitrariness of the initial flashes. Of course, the mathematics of 
the model works with any choice of the initial flashes. 

3.4 Hyperboloids and Cauchy Surfaces 

A hyperboloid need not be a Cauchy surface in ^i] in fact, it never is in Minkowski 
space-time,^ though it is "almost Cauchy" in the sense that most complete timelike 
curves do intersect a given hyperboloid: the ones avoiding the hyperboloid accelerate to 
the speed of light. However, it is not necessary for our purposes that the hyperboloids 
be Cauchy surfaces. It is sufficient that, for a single particle, 

the Dirac equation defines a unitary evolution operator 

Uf : L2(E) ^ L2(^) for every Cauchy surface E (25) 

and every hyperboloid Jif. 



^To see this, consider the following example of a complete timelike curve x{t) that does not intersect 
the unit hyperboloid J^(0, 0,0, 0): the uniformly accelerated curve x{t) = (ci,0, 0,-\/l + c^f^)- 



We conjecture, but do not have a proof, that this is the case in Minkowski space-time 
for a large class of vector potentials A^. We hope to be able to provide a proof in a 
future work. Here is a partial argument: To show that the evolution is unitary requires 
to show that it is onto and norm-preserving. To show the latter it should suffice to 
show that, for all wavefunctions ijj from a dense subspace of L^(S), no more than a set 
of measure zero of flow lines of the IV'P distribution (Bohmian trajectories) avoid the 
hyperboloid. This presumably amounts to the vanishing of the probability flux into the 
future null infinity. This could be expected to be the case for ip from the form domain 
of the Hamiltonian, because to accelerate a positive amount of Bohmian trajectories to 
the speed of light should require, in some sense, an infinite amount of energy, while such 
ip has finite energy expectation. Indeed this conclusion apparently follows, for static 
electromagnetic fields tending to zero fast enough at infinity and for a suitable class of 
t/j's, from the fiux-across-surfaces theorem J3 using the global existence of Bohmian 
trajectories [2H]- 

How about other space-times than Minkowski? In some space-times, hyperboloids ac- 
tually are Cauchy surfaces and so the problem is absent; for example, think of Minkowski 
space-time modulo a spacelike 3-lattice. Relevant conditions for this case may be empti- 
ness of the future null infinity and absence of future singularities. The future null infinity 
could be expected to be empty when "space has finite volume growing not too quickly," 
such as in Minkowski space modulo a spacelike 3-lattice. In space-times in which, as in 
Minkowski space-time, hyperboloids are not Cauchy surfaces, a relevant condition for 
fl25|) seems to be that, for all ip from a dense subspace, the probability fiux into the 
future null infinity vanishes. 

For particles with zero rest mass, such as photons and gravitons, one would expect 
that there typically is a positive probability fiux into the future null infinity.^ One 
way of avoiding this problem would be to postulate that fiashes are associated only 
with massive particles, and to have the arguments of the wavefunction corresponding to 
massless particles get integrated out in ^lU^ along an arbitrary Cauchy surface. 

3.5 Nonlocality 

The model is nonlocal. This is manifest in ()19|) in that the joint distribution of the fiashes 
Yi and Yj, j ^ i, does not factorize. One can easily find situations in which the events 
Yi and Yj (now taken to lie in the same space-time manifold) are spacelike separated, 
and still the distribution of Yi depends on the realization of Yj] or, of course, we may 
view it the other way round: that the distribution of Yj depends on the realization of 

3.6 Generations 

Since the fiashes were constructed here together in generations, i.e., groups Yi, . . . ,Yn oi 
N fiashes, one may wonder whether the theory presupposes, or provides, some structure 



^This remark has been kindly pointed out to me by Fay Dowker. 



beyond the mere flashes and their labels, a structure defining which fiashes belong to the 
same generation. Such a structure would seem against the spirit of relativity, as much 
perhaps as a preferred slicing of space-time into hypersurfaces, since it would define, 
for two spacelike separated fiashes, which is earlier and which is later, if they belong to 
different generations. 

In fact, however, the theory neither presupposes nor provides a grouping of the flashes 
into generations. Rather, one could take as the initial flashes a combination of some of 
the flrst-generation flashes Xj and some of the second-generation flashes Yi, and obtain 
the same (conditional) distribution of the future flashes. This will become manifest in 
the reformulation of the law of the flashes of Section 14.11 

3.7 Negative Energy Contributions 

In the context of the Dirac equation, one often considers wavefunctions with negative 
energy as unphysical, more precisely those wavefunctions containing, for at least one 
particle, at least some contribution from the negative-energy subspace of the 1-particle 
Dirac Hamiltonian. Such wavefunctions we will call "non-positive" for short, and the 
other wavefunctions, containing exclusively positive-energy contributions, "positive." 
The collapse accompanying a flash will generically map a positive wavefunction to a 
non-positive one. This fact poses a problem if we want to take the model seriously, 
and in particular if we want to extend it to quantum fleld theories incorporating anti- 
particles. 

One might conclude that the law for collapsing the wave function should be so mod- 
ifled that the collapsed wavefunction remains positive, for example by simply projecting 
the wavefunction, on top of the usual multiplication by a Gaussian as in ()23|) . to the 
space of positive wavefunctions. In formulas, if P+ denotes this projection, and if we 
abbreviate the collapse law (j^ as 

$ oc j\[^ , 

the modifled collapse law corresponds to 

$ oc P+iJ"^) . (26) 

Since P+ is the same operator as convolution with a suitable function whose width is 
small (of the order of magnitude of the Compton wavelength), and thus does not change 
the |$p distribution a lot, replacement of (j^^ with (f^ would perhaps entail merely a 
small change in the distribution of the flashes. 

Alternatively, one might conclude that the spontaneous collapses can lead to spon- 
taneous pair creation, since non-positive wavefunctions would appear to have something 
to do with anti-particle states. This effect might entail observable deviations from quan- 
tum mechanics. For calculating concrete predictions, it would seem necessary to flrst 
precisely formulate a version of the model that incorporates particle creation and uses 
wavefunctions from Fock space, but one might guess already from the present model 
that the probability of pair production at a collapse could be of the order of magnitude 

of||(l-P+)(j^)||Vllj^f- 

10 



4 Distribution Formulas 

In the previous section we have defined, in a rather abstract way, the distribution of 
the fiashes in space-time. We will find it helpful to have further expressions for this 
distribution; such expressions we provide in this section. 

4.1 The First n Flashes 

It is now straightforward, though somewhat tedious, to derive an explicit formula for the 
joint distribution of the first rii fiashes of type i, given the "initial" fiashes X°, . . . , X^. 
For any choice of spacelike (Cauchy) surfaces S^, we can express the distribution in 
terms of the initial wavefunction \E'°, restricted to the Cartesian product of the S°, in 
the form 



Probfe G rf^x(=; i = 1, . . . ,N-k = 1, . . . ,n^ = /"^^ E^") ( JJ JJ d^x^ J 



vl/0), (27) 



^i=i fc=i 



where the scalar product is taken in (8)jL^(Il°). Here, n = {rii, . . . ,nN), and E^"^^ is a 
positive-operator- valued measure (POVM) on Yli-^i^'- ^^ i^ of a product form. 



N Hi s N 






^i=\ k=\ ^ 1=1 ^fc=l 

where Ef^ = E^^^J is a POVM on ^", defined on L^T.^^) by 

^A:=l ^ ^fc=l ^ 

' ' ^ k=i ^ 



(29) 



X T^ exp ( -— ;> ^ t-clist(x;% x'l ^ j ) 3\ f, ■ ■ ■ j7 j7 ■ ■ ■ f, j] 



with 1b the indicator function of the set B and j^ the self-adjoint "collapse" operator 
defined on L^(S°) as follows. Let, for any spacelike surface S and a; G S, js(a^) be the 
multiphcation operator on i^^(E) that multiplies by the function Jt;{x, •) defined in (|7I|). 
Then 

In our notation E^J we conceal the (uninteresting) dependence on the choice of S°. 
Normalization of -E^ , i.e., E'j (^") = 1, can also directly be seen from the fact that 

Et'-'U\\d^x^xnx7)]=Er{\\d'x^], (31) 



f[d'x'lxnx-)\=El-^ff[d'x'l), 

k=l ' ^k=l ' 



11 



which in turn follows from 



J?(x) 



d^y—exp{^ t-dist(|/,a;) ) jji^'(y,^)iyy 



1. 



(32) 



Eq. (j77|) is, in effect, a Heisenberg picture formulation, as illustrated particularly by 
the way the unitary operators occur in (jHUj) . and by the arbitrariness of the surfaces S°, 
which may even lie in the future of (some of) the flashes. 

4.2 The Flashes up to Given Surfaces 

It is a corollary of (fTTj) that the probability for obtaining, up to spacelike surfaces Sj, a 
particular sequence xj, . . . , x"' of flashes for every type i, is 



Probfe n ^(Si) G {x°} X d^x] X ... X d'^x^ 



x?^ Vi 



¥ 



¥ 



(33) 



£;(-+!) (JJd^xl X ■ ■ ■ X d^xl' X ^(E,) 

i 

with n + 1 := (ni + 1, . . . , tzat + 1). 

4.3 Flash Rate in the Temporal Picture 

It is often desirable to select a time coordinate — a slicing parametrized by t into spacelike 
surfaces Ej(t) — in each ^i and to employ a picture in which everything evolves as a 
function of the time coordinate. In this picture, the interesting quantity is the probability 
of having a flash of type^ / G {!,..., A^} between t and t + dt, conditional on the flashes 
up to time t. For this we obtain the following formula, in which we assume that the 
given flashes are timelike separated, x^ G ^(xf~^) for all i and k < rii, and that d^y is 
a volume element in S/(t). We denote by dt x d^y the 4- volume element between S/(t) 
and ^/(t + dt) swept out by the normals on S/(t) over d^y; the volume of this element 
is c{y) dtd^y with c{y) dt = t-dist(S/(t + dt),y). 



FrohfOj n ^(S/(t)) n ^(S/(t + dt)) = {Y}, Y edtx d^y 

Q,n^{^,{t)) = {xl...,xr}yt) = 

= c{y)dtd^y l,^r^"i)iy) — exp( t-dist(|/,x"^ 



(34) 



X 



^i 



^2^(^(^iW))®---®^"^(^)'®--- 



^(1) 



^;>(^(s-w) 



^. 



^y 



(8),Efj„,(^(s.(t)) 



^y 



^Unlike as usual, the capitalization is not meant here to indicate that / is a random variable, but 
merely to distinguish that particular number from the other z's. 
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where 

^t = lt(6djT---jlW (35) 



with (arbitrary^) normahzation factor 7^, and ji{y) (appearing in the /-th factor of the 
tensor product) is the collapse operator corresponding to a flash at y, 

5 The Temporal Picture 

We now describe the model in the temporal picture, based on a slicing of space-time. 
This is basically a discussion of (jHH), the flash rate formula of the temporal picture. 
With every time t we associate a wavefunction 

** = 7*((8)f>g^*^jT---jl)^°, (37) 



the same as in (jHSj) except that we now shift from the Heisenberg to the Schrodinger 
picture and use the unitary Dirac evolution to define ^t on Si(t) x ■ ■ ■ x SAr(t). ^t 
evolves unitarily apart from collapses. We remark that ^t in fact depends only on the 
surface at time t: it depends on the slicing of space-time only through this surface, in 
the sense that different slicings that happen to coincide at t have the same ^f As is 
clear from (jH^j) . ^t is determined by the initial wavefunction '^'^ and the flashes up to 
time t. 

Whenever a flash occurs, the wavefunction \&j gets collapsed by applying to it the 
suitable collapse operator j"'^ (and then renormalizing) . Recall that also in the GRW 
model, the wavefunction gets collapsed by applying the appropriate collapse operator. 
However, while in the GRW model the collapse operator is always a multiplication oper- 
ator in the position representation, this is not necessarily the case in our model; instead, 
j"'^^ arises from a multiphcation operator on the suitable hyperboloid, ^(X"*"*"^, X"*), 
by using the Dirac propagators to get to L^(^(Xf''^ i^T)) ^"^^ back; in other words, 
j"'"^ is a multiplication operator evolved, as in the Heisenberg picture, to another sur- 
face. 

To compute the flash rate, it is unnecessary to know the entire history of collapses. 
Instead, it suffices to know the present wavefunction ^^ and the last collapses of all 
types, Xf\ In other words, the evolution of the (A^ + l)-tuple (^j,X"\ . . . ,X]^^) is 
a Markov process. Although this is true as well in the GRW model, there, by way of 
contrast, already '^t itself follows a Markov process, and knowledge about the space-time 
locations of the previous flashes does not add any information, beyond what is encoded 
in \Eft, about the distribution of the future flashes. 

^That is, although we will consider the usual normalization {'^t\'^t) = 1, l|34|l holds for any choice 
of 7t, for instance 74 = 1. Another useful possibility is to choose 7^ such that the denominator in the 
last line of H34(l becomes 1. 
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Another difference between the ffash rate formula of our model ()34|) and that 
of the GRW model (jH)) is the appearance of a denominator and an further factors, 
£^j- (^(Ej(t))), in the numerator; they are needed to correct the rates for conditioning 
on the knowledge that no further collapses have occured until time t. 

How does the wavefunction transform under a change of slicing of space-time, such 
as given by a Lorentz boost in Minkowski space-time? In two ways. First, some flashes 
may lie in the future of the new surface S^(t) but in the past of the old surface Sj(t) and 
vice versa; consequently, the corresponding wavefunctions differ by application of the 
collapse operators (respectively their inverses, and renormalization) belonging to these 
flashes. Second, on top of that the wavefunctions differ by the unitary Dirac propagator 
from one surface to the other. 

6 The Low Velocity Regime 

In this section we consider Minkowski space-time and assume the existence of a Lorentz 
frame in which all velocities, at which wave packets move or spread, are small compared 
to the speed of light. We point out that in this regime, and in this frame, our model 
approaches the GRW model. 

The basic observation is that for large t, such as t of the order of magnitude of 
r ^ 10^^ sec, the hyperboloid ^ct becomes rather fiat and is well approximated by 
the plane x° = ct tangent to the hyperboloid, at least in the neighborhood of their 
common point (ct, 0, 0, 0) with radius R such that R <^ ct, or R <^ 10*^ light-years. As a 
consequence, to multiply by a Gaussian with width a <^ R centered somewhere inside 
this neighborhood means approximately the same on the hyperboloid and on the plane. 
And since, by the assumption on velocities, the wavefunction will lie almost completely 
within this neighborhood for t of the order of magnitude of r, the action of the collapse 
operator agrees approximately with that of the GRW model. Similarly, the factor l,^(x) 
is 1 on the support of \E'j. 

Replacing our collapse operators with those of GRW, we obtain E^ ■^(^(Sj(t))) ^ 
^-{t-Ti)/T 1^ gQ that, since c{y) = c for this shcing, (jM|) in fact reduces to (jHI). 

There is a subtlety in that if N is very large, flashes do sometimes occur for which 
the timelike distance cAT from the previous flash of the same type is not of the order 
of magnitude of ct, but much smaller. (However, cAT is still much larger than a except 
in very rare cases that occur only once in cr/a ~ 10^° collapses.) Fortunately, relevant 
wave functions obey a bound on their spread: it is at most a initially (thanks to the 
previous collapse) and grows much slower than at rate c. Thus, the spread is still 
much smaller than cAT, so that the collapse operators of the two models do not differ 
appreciably on these wavefunctions. 
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7 Predictions 

What predictions does the model entail? To what extent are they in agreement with 
quantum mechanics? In what way do they deviate from the predictions of the nonrela- 
tivistic GRW model? 

To begin with, a difficulty with obtaining any predictions at all from the model is 
that it does not involve any interaction, and thus does not support the formation of 
macroscopic bodies such as observers or apparatuses. However, we can say what the 
predictions will be like once interaction is included in whatever way, be it by particle cre- 
ation and annihilation or by a modification of the unitary propagators: we can consider 
a wavefunction as would arise from interaction. 

As a corollary of the previous section, for any experiment in which no parts need 
move at relativistic speeds our model approximates the GRW model and is thus in 
agreement with quantum mechanics to the same extent as the GRW model; this in- 
cludes all presently doable experiments; for a discussion of future experiments that may 
distinguish spontaneous collapse theories from quantum mechanics, see Section V of jHj- 
The model is in particular in agreement with the result of EPR-Bell experiments such 
as Aspect's |2j and thus violates Bell's locality inequality |H]. 

A general pattern of behavior of the model follows from another trait it has in 
common with the GRW model: that disentangled subsystems are governed by the same 
laws as the whole, and follow an independent collapse process. Therefore, a small 
system of A^i < 10^ particles will not collapse for the next thousand years, provided it 
stays disentangled. A macroscopic body, however, cannot support superpositions over 
distances much wider than a ^ 10^^ m for longer than a split second. 

The most obvious deviation of our model from the GRW model is that a system 
moving at a speed close to c will have a reduced rate of spontaneous collapse, reduced 
by just the factor that one would expect from a naive application of time dilation. 
This will be hard to see in experiment, of course, given that as yet we cannot see any 
spontaneous collapses in experiment. 

Extensions of the model to quantum field theory might deviate from the nonrelativis- 
tic GRW model in a prediction of spontaneous pair creation at the fiashes, as discussed 
in Section 1X71 

Superluminal communication on the basis of entangled sets of particles is impossible 
in our model. To see this, suppose that Alice and Bob are widely separated and share 
a system of entangled particles; the marginal distribution of the Qi for the z's of all 
particles located on Alice's side is, as follows from (fTfj) and (j^ . independent of the 
fields applied to the Xi for the i's of the particles on Bob's side, such as the metric and 
electromagnetic vector potential of ^j. 

8 Perspective 

Our model seems to be the first model in the literature that achieves all of the following: 
(i) it describes a possible (many-particle) world in which outcomes of experiments 
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have (to a sufficient degree of accuracy for all cases presently testable) the proba- 
bilities prescribed by quantum theory, 

(ii) it describes objective events in space-time (the flashes), in contrast to theories 
merely associating a wavefunction with every spacelike surface, 

(iii) it is fully compatible with relativity in that it does not rely on a preferred slicing 
of space-time, and 

(iv) it works in the continuum, in contrast to theories assuming a discrete space-time. 

We give a brief overview of the literature concerning relativistic models explaining 
the probabilities prescribed by quantum theory. Such models come in two varieties, 
either as a variant of Bohmian mechanics or as a spontaneous collapse theory. 

Among the variants of Bohmian mechanics, Bohm jH] gives a Lorentz-invariant equa- 
tion of motion for a single Dirac particle; Bohm and Hiley ^ give a many-particle version 
based on a preferred Lorentz frame; Diirr et al. |16| generalize to an arbitrary spacelike 
preferred slicing of space-time, possibly determined by a covariant law involving the 
wavefunction; Samols [23 I2Z] gives a Bohm-type model on a discrete space-time using 
a preferred slicing; Goldstein and Tumulka [201 give a nonlocal many-particle version 
without preferred slicing, which however fails to yield any probabilities. Both Berndl et 
al. [7j and Dewdney and Horton 10^ suggest, instead of a preferred slicing, a preferred 
joint parametrization (or synchronization) of the world lines, and thus obtain a nonlocal 
Bohm-type dynamics; however, this does not really conform any better with the spirit of 
relativity than a preferred slicing; in addition, the models fail to yield any probabilities. 

All spontaneous collapse models deviate slightly from the quantum prescriptions 
in their probabilities. An overview of spontaneous collapse models is given by Bassi 
and Ghirardi in [3]. Dove and Squires [121 US] have made steps towards a relativistic 
model based on discrete flashes.^ Dowker and coworkers ^5 |^ give a collapse model 
on a discrete space-time that does not need a preferred slicing; it is not known how 
this model could be adapted to a continuum. All other efforts towards a relativistic 
collapse model are based on the approach of continuous spontaneous localization (CSL), 
corresponding to diffusion processes in Hilbert space; some references describing research 
in this direction are |21 CHI IH 123 123 • 

A somewhat surprising feature of the present situation is that we seem to arrive 
at the following alternative: Bohmian mechanics shows that one can explain quantum 
mechanics, exactly and completely, if one is willing to pay the price of using a preferred 
slicing of space-time; our model suggests that one should be able to avoid a preferred 
slicing if one is willing to pay the price of a certain deviation from quantum mechanics. 



^Their paper ll^l (reprinted, with minor extensions, as Chapter 7 of ,12; ) is irritating in that the 
authors claim to provide a Lorentz-invariant collapse model but do not keep their promise. They define 
an evolution of the wavefunction given the flashes, according to which the wavefunction collapses on 
the future light cone of each flash and evolves unitarily in between. However, they do not specify a 
probability law for the flashes. 
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Note added. Several articles discussing the relativistic model of this paper have been 
written since its first preprint version became available in 2004 [^ 1221 P . 
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